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We discuss the role of the conservation laws related with U(l) internal symmetry group in the 
statistical model description of particle productions in ultrarelativistic heavy ion collisions. We 

, derive and show the differences in particle multiplicities in the canonical and the grand canonical 

formulation of quantum number conservation. The time evolution and the approach to chemical 
equilibrium in the above ensembles is discussed in terms of kinetic master equation. The application 
of the statistical model to the description of (multi)strange particle yields at GSI/SIS and the SPS 

' energies is also presented. 

o : 

\D ! I. INTRODUCTION 

. There are in general two approaches to describe integrated particle yields measured in ultrarelativistic heavy ion 
collisions: (i) the microscopic transport models and (ii) the macroscopic statistical thermal models [|]-|]]. ^ n * n ' s 
f — , article we will discuss the statistical approach and show that it provides a very satisfactory description of experimental 
data. We will emphasize the importance of the conservation laws in the particular strangeness conservation when 
CO modelling particle chemical freczcout conditions. 

Within a statistical approach, the production of particles is commonly described using the grand canonical (GC) 
ensemble, where the charge conservation is controlled by the related chemical potential. In this description a net 
value of a given U(l) charge is conserved on the average. The (GC) approach can be only valid if the total number 
of particles carrying quantum number related with this symmetry is very large. In the opposite limit of a small 
particle multiplicities, conservation laws must be implemented exactly and locally, i.e., the canonical (C) ensemble 
for conservation laws must be used The local conservation of quantum numbers in the canonical approach 

severely reduces the phase space available for particle productions. This treatment of charge conservation is of crucial 
importance in the description of particle multiplicities in proton induced processes [p|jic|l, in e + e || as well as in 
central heavy ion collisions at low beam energies Oj . 

In this article we describe the exact strangeness conservation in the context of relativistic statistical thermodynamics. 
A kinetic theory for the time evolution of particle production and the approach to the grand canonical and the canonical 
equilibrium distribution will be also introduced. Finally the example of the applications of the statistical model in 
■ (C) ensemble is presented in the context of low energy central as well as in high energy peripheral heavy ion collisions. 
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II. STATISTICAL MODEL AND PARTICLE MULTIPLICITY 



The exact treatment of quantum numbers in statistical mechanics has been well established for some time now 
JfyplJIl]]. ^ i s m general obtained by projecting the partition function onto the desired values of the conserved charges 
by using group theoretical methods. For our purpose we shall only consider the conservation laws related to the 
abelian U(l) symmetry group. In particular, we concentrate on strangeness conservation. 

The basic quantity in the statistical mechanics describing a thermal properties of a system is the partition function 
Z(T,V). In the (GC) ensemble, 

Z GC (T, V, fig) = Tr[e-^ H ~^] (1) 

where Q is the conserved charge, H the hamiltonian of the system, \iq is the chemical potential which plays the role of 
the Lagrange multiplier which guarantees that the charge Q is conserved on the average in the whole system. Finally 
(3=1 /T is the inverse temperature. 
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In the (C) ensemble the charge Q is conserved exactly. Thus there is no more chemical potential under the trace 
and instead we calculate the partition function summing only these states which are carrying exactly the quantum 
number Q, that is 

Z^(T,V)=Tr Q [e^ a ] (2) 

The canonical and the grand canonical partition functions are related through the following cluster decomposition, 

+00 

Z GC {T,V,\)= J2 xQz q( t ,V) (3) 

Q=— 00 

where the fugacity A = exp(/3//g) and the sum is taken over all possible values of the charge Q. 

For the (GC) partition function, which is well behaving analytic function of the fugacity A, the above relation can 
be inverted and the canonical partition function with a given value of the charge Q reads, 



1 
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X: 



-iQ<p 



Z{T,V,J>) (4) 



where the generating function Z is obtained from the grand canonical partition function replacing the fugacity 
parameter A by the factor e l< ^, 

Z(T,V,(f>) = Z GC {T,V,\^ e^) (5) 

The form of the generating function Z in the above equation is model dependent. Having in mind the application of 
the statistical description to particle production in heavy ion collisions we calculate Z in the ideal gas approximation, 
however, including all particles and resonances H]. This is not an essential restriction, because, describing the freeze- 
out conditions we are dealing with a dilute system where the interactions should not influence particle production 
anymore. We neglect any medium effects on particle properties. In general, however, already in the low-density 
limit, the modifications of resonance width or particle dispersion relations could be of importance |lS^| . For the sake 
of simplicity, we use classical statistics, i.e. we assume temperature and density regime so that all particles can be 
treated using Boltzmann distributions. 

In nucleus- nucleus collisions the absolute values of the baryon number, electric charge and strangeness are fixed 
by the initial conditions. Modelling particle production in statistical thermodynamics would in general require the 
canonical formulation of all these quantum numbers. From the previous analysis ffljll||, however, it is clear, that 
in heavy ion collisions only strangeness should be treated exactly, whereas the conservation of baryon and electric 
charges can be described by the appropriate chemical potentials in the grand canonical ensemble. 

Within the approximations described above and neglecting first the contributions from multi-strange baryons, the 
generating function in equation (|j) has the following form, 

Z(T, V, fx Q , [is-, 4>) = exp(7V s=0 + N s=1 e^ + N^-ie'**) (6) 
where iV s= o,±i is defined as the sum over all particles and resonances having strangeness 0,±1, 

k 

and Z\ is the one-particle partition function defined as 

Zl = ^^mlTK 2 (m k /T) exp(b k ^ B + qkfJ-q) (8) 

with the mass rrik, spin-isospin degeneracy factor g k , particle baryon number bk and electric charge q k . The volume 
of the system is V and the chemical potentials related with the electric charge and the baryon number are determined 
by \i q and fis respectively. 

With the particular form of the generating function equations (^|,0J^) the (/^integration in equation (|4|) can be done 
analytically giving the canonical partition function for a gas with total strangeness S [[ll] : 

Zs{T, V, hb, Mq) = z o(T, V, ^ b ,^q)Is{x) (9) 
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where Zq = exp (Ns=o) is the partition function of all particles having zero strangeness and the argument of the 
Bessel function 



x = 2y/S 1 S- 1 . (10) 

with S±i = N s= ±i. The parameter x thus measures the total number of strange particles in thermal fireball. 

The calculation of the particle density in the canonical formulation is straightforward. It amounts to the 
replacement 

Z\ i * x k zl (11) 

of the corresponding one-particle partition function in equation (Q) and taking the derivative of the canonical partition 
function equation (|l|) with respect to A& 

n^ = [Xk-^lnZ Q (X k )] Xk=1 (12) 

As an example, we quote the result for the density of thermal kaons in the canonical formulation assuming that the 
total strangeness of the system S = 0, 



Z\ Si h(x) 
V y/s7S~[ Iq{x) 



(13) 



Comparing the above formula with the result for thermal kaons density in the grand canonical ensemble, = 
(Z^/V) exp (/i s /T), one can see that the canonical result can be obtained from the grand canonical one replacing the 
strangeness fugacity Ag = exp (fig/T) in the following way: 



C _ GC 
l K — n K 



Si h(x) 



(14) 



In the limit of large x that is large volume and/or temperature the canonical and the grand canonical formulation 
are equivalent. For a small number of strange particles in a system, however, the differences are large. This can be 
seen in the most transparent way when comparing two limiting situations: the large and small x limit of the above 
equation. In the limit noowe have 

x— »oo 1q(X) 

and the kaon density is independent of the volume of the system as expected in the grand canonical ensemble. On 
the other hand in the limit of a small x we have 

km - " (16) 

and the particle density is linearly dependent on the volume. It is thus clear, that the major difference between the 
canonical and the grand canonical treatment of the conservation laws appears through different volume dependence 
of strange particle densities as well as strong suppression of thermal particle phase space. The relevant parameter, 
Fs, which measures the suppression of particle multiplicities from their grand canonical result is determined by the 
ratio of the Bessel functions 

Fs = (17) 

with the argument x defined in equation (|l0|). 

In Fig. 1 we show the canonical suppression factor Fs(x) as a function of the argument x. To relate the initial volume 
of the system with the number of participant in A-A collisions one uses the approximate relation V ~ l.9TrA part . 
The corresponding value of x at SIS, AGS and SPS energy is calculated with the baryon chemical potential and 
temperature extracted from the measured particle multiplicity ratios p"l| ] . The results in Fig. 1 shows the importance 
of the canonical treatment of strangeness conservation at SIS energy. Here, the canonical suppression factor can be 
even larger than an order of magnitude. 
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For central Au-Au collisions at AGS or SPS energy this suppression is not relevant any more and the (GC)-formalism 
is adequate. In general, one expects, that the statistical interpretation of particle production in central heavy ion 
collisions requires the canonical treatment of strangeness conservation if the CMS collation energy y/s < 2— 3GeV. This 
is mainly because at these energies the freeze-out temperature is still too low to maintain large argument expansion 
of the Bessel functions in equation ([l5|). The canonical description of strangeness conservation can be, however, 
also of importance at the SPS energy were ^/s ~ 18GeV when one considers the peripheral heavy ion collisions. 
This is particularly true for multistrange particle production since the canonical suppression of the thermal particle 
phase-space increases with strangeness content of the particle pi] . 
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FIG. 1. Canonical strangeness suppression factor (see text). 



A. Multistrange particle multiplicities 

The extension of the canonical description to multi-strange particle multiplicities is straightforward. One needs 
first to extend the generating functional in Eq. 6 by including the contributions of multistrange baryons. In this case 
the canonical partition function constraint by the strangeness neutrality condition reads [|lOUl3| , 



Zs=o = ^- # exp(]T S s e^), 



(18) 



where S s = J^i %l an d the sum is taken over all particles and resonances carrying strangeness s. The one-particle 
partition function Z\ is defined in Eq.(8). 

With the particular form of the partition function given by Eq. (|l8|) the density n s of particle i with strangeness s 
in volume V is obtained by the replacement Zi t— > \Zi in Eq. ( |l8| ) and then taking an appropriate derivative ]7|^0|]: 

(N± a ) AdlnZn (S ±1 ) s 

n ±s = T , = [tt^t; — jAi=i - Z± 



oo 

{/.(a:i)/o(a:2) + ^ Irn{x 2 )[I 2mM 4 Xl )A m l 2 + I 2m ^ s] { Xl )A' m ' 2 ]} /Z S=0 (19) 



where 



sfsl 



x k = 2^ST k , A=-^, (20) 



and the partition function 
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oo 

Z s =o * Io(xi)I (x 2 ) + h m {xx)I m {x2)[A m ' 2 + A'™/ 2 }. (21) 

m— 1 

In the derivation of Eg. (|l9|) we have neglected, after differentiation over particle fugacity, the term S±s. This 
approximation, however, due to small value of S±3 coefficients in comparison with S±± and S±2 is quite satisfactory. 
The complete expression for the partition function without any approximation can be found in reference p3| . 

In the large system like Pb-Pb and for large collision energy, required to reach high T, the density n s of particle 
carrying strangeness s is V independent. In the opposite limit, however, this dependence is changed to n s ~ V s , 
which can be verified from Eq. (|T^). Indeed for small Xi we have approximately: 

n± s ~ Z ±s — — , — — r. (22) 
(5 + i5_i) s /^ Io(xi) 



and when expanding the Bessel functions I n (x) ~ x n in Eq. (|22j) we see that n s ~ V s . 

From the above expression it is clear that strangeness suppression, which is measured by the ration I s /Iq, is 
increasing with the strangeness content of the particle. Thus, there are two important ingredients of canonical 
modifications of multistrange particle density with respect to their (GC) value: (i) the density is volume dependent that 
is also centrality dependent and (ii) the thermal phase is suppressed and this suppression increases with strangeness 
content of the particle. 




FIG. 2. Measured K + multiplicity per participant A par t as a function of Apart for Au-Au collisions at 1 A/GeV 
together with the canonical model results 

III. TIME EVOLUTION AND STRANGENESS EQUILIBRATION 



In the last section we have formulated the statistical model for strange particle multiplicities < Ns > assuming that 
the system is in thermal and chemical equilibrium. We have shown that dependently on the total number of strange 
particles there are two distinct equilibrium limits |qjll]|: if < N s > is small then we are in the canonical regime and 
in the opposite limit the canonical and grand canonical description coincides. In this section we consider the time 
evolution of the multiplicity of strange particles and formulate a kinetic master equations which distinguish between 
these two equilibrium limits. For a sake of illustration we consider a simple example of K + K~ production in the 
environment of thermal pions in volume V and temperature T due to the following binary process, 7r + 7r~ — * K + K~ . 
We formulate for this example the kinetics for the time evolution of kaon multiplicities and their approach to chemical 
equilibrium . 

In the standard formulation Jl3] , the rate equation for this binary process is described by the following population 
equation: 



d < N K > G Ar 



>< N. 



K~ 



>, 



(23) 
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where G = (ctgv) and L = (<jlv) give the momentum-averaged cross sections for the gain 7r + 7r~ — > K + K~ and the 
loss K + K~ — > 7r + 7r~ process respectively. The value of < Nk > represents the total number of produced kaons. 

To include the possible correlations between the production of K + and K~ B, let us define Pij as the probability 
to find i number of K + and j number of K~ in an event. We also denote by Pj as the probability to find i number 
of K pairs in an event. The average number of K per event is defined as: 

oo 

(N K )=J2iPv (24) 

i=0 

We can now write the following general rate equation for the average kaon multiplicities: 

i,j 

Due to the local conservation of quantum numbers, we have: 

p. . — p x. . 

X'.' 7 ' - = E 1 ^ = ^ = ( jV ) 2 + (^ 2 )' ( 26 ) 
»,j » 

where (SN 2 ) represents the event-by-event fluctuation of the number of K + K~ pairs. Note that we always consider 
abundant 7r + and 7r~ so that we can neglect the number fluctuation of these particles and the change of their 
multiplicities due to the considered processes. 

Following Eqs.(|2^-p6[) the general rate equation for the average number of K + K~ pairs can be written as: 

^ = fWrOWr-) - (27) 

For abundant production of K + K~ pairs where (Nk) ^> 1, 

(N 2 K ) « (iV^) 2 , (28) 
and Eq. ( p7| ) obviously reduces to the standard form: 

^«f<^>(JV>-^W. (29) 



However, for rare production of K + K pairs where (Nk) <1, the rate equations ( P3[ ) and (29) are no longer valid. 
We have instead 

(N 2 k ) « (JVjc), (30) 

which reduces Eq.(p7|) to the following form 

^ « ^<^><^-> - ^W- (3D 

Thus, in the limit where (Nk) *C 1, the absorption term depends on the pair number only linearly, instead of 
quadratically for the limit of (Nk) 3> 1. Thus, it is clear that the time evolutions and equilibrium values for kaon 
multiplicities are obviously different in the above limiting situations. 

In the limit of large < Nk > the equilibrium value for the number of K + K~ pairs, which coincides with the 
multiplicity of K + and K~ , is obtained from Eq.(29) as, 

< Nk > G q C = ^m K TK 2 (M K /T) (32) 

thus, it is described by the (GC) result with vanishing chemical potential due to strangeness neutrality condition. 

In the opposite limit where (Nk) <C 1, the time evolution is described by Eq.(|3l]), which has the following equilibrium 
solution: 
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N c = 

1 Y cq 



^M 2 K+ TK 2 (M K+ /T) 



X_M 2 K -TK 2 {M K -/T) 



(33) 



The above equation demonstrates the locality of strangeness conservation. With each K + the K~ is produced 
in the same event in order to conserve strangeness exactly and locally. This is the result expected from the (C) 
formulation of the conservation laws as described in the previous section. We note that Eq.(|33|) is just the leading 
term in the expansion of the canonical result for multiplicities of particles which are carrying U(l) charges. The 
general expression is given by Eq.(13) and Eq.(20). 

Comparing Eq.(32) and Eq.(33), we first find that, for (Nk) <C 1, the equilibrium value in the canonical formulation 
is far smaller than what is expected from the grand canonical result as 



<N K > 



eq 



\<N K > 



GCi 
oq J 



< < N K > 



GC 
eq ■ 



(34) 



This shows the importance of the canonical description of quantum number conservation when the multiplicity of 
particles carrying non-zero U(l) charges is small. We also note that the volume dependence in the two cases differs. 
The particle density in the GC limit is independent of V whereas in the opposite canonical limit the density scales 
linearly with V. Secondly, we note that the relaxation time for a canonical system is far shorter than what is expected 
from the grand canonical result ||. It is also clear from Eq. (|26|) that (C) and (GC) limits are essentially determined 
by the size of (8N 2 ), the event-by-event fluctuation of the number of K + K~ pairs. The grand canonical results 
correspond to small fluctuations, i.e., (SN 2 ) / (N) 2 « 1; while the canonical description is necessary in the opposite 
limit. 



IV. MODEL PREDICTIONS VERSUS EXPERIMENTAL DATA 



In heavy ion collisions the number of produced strange particles depends on the collision energy and the centrality 
of these collisions. At low collision energies like eg. in GSI/SIS the freezeout temperature is relatively low being of 
the order of 50 — 80MeV. Consequently the number of strange particles in the final state is very small. Following our 
discussion in the last sections it is clear that the statistical description require here the canonical approach. In Fig. 2 
we show the experimental data on K + yield per participant A part as a function of A part measured in Au-Au collisions 
at Ei ab — 1 A/GeV @. 



\ ' ' ' ' i 

Ni -Ni 7.9 A GeV 



1 



,f?=4 . 2fm 
sR = 4 . Ofm 

,-' / 
.' y --R=3 8fm 

.'' / ,'' 
,-' / ,'' 

,'' y --' 
■■ / .-' 
x y 
1 



/' 

,-' / 



70 75 
T[GeV] 



FIG. 3. The ratio of S to K + multiplicity as a function of temperature. The point with errors indicates the predictions 
of the thermal model for Ni-Ni collisions at 1.9 A/GeV. The lines corresponds to canonical calculations for different 
centrality measured by the initial radius R of the system. 

The volume parameter in the statistical model scale with the number of participant. We can thus directly compare 
the model with data by fixing thermal parameters, the temperature and the baryon chemical potential, such that 
to reproduce the measured particle multiplicity ratios In Fig. 2 the results of the canonical model is shown 
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by the full line. The results clearly indicates that strong almost quadratic dependence of kaon yield on the number 
of participant is well reproduced by the model. The quadratic dependence of particle multiplicity in the canonical 
regime is the basic property of the model. In Fig. 3 we calculate the ratio E/K + multiplicities for central Ni — Ni 
collisions in the temperature range which corresponds to the collision energy of Ei a b ~ 2 A/GeV. 

The yield of 2 is seen to be substantially smaller then the yield of K + . This result is not only related with the 
differences in particle masses but particularly it appears due to the canonical suppression. Since S carry strangeness 
minus two it has to be produced together with two K + to locally and exactly neutralized strangeness. It is also 
interesting to note that the S/_fT + ratio is independent on the baryon chemical potential. This is because K + appears 
together with A, thus it contains the same dependence on as multistrange baryons. 

The importance of the canonical treatment of strangeness conservation is also seen in higher collision energies 
like at the SPS when considering centrality dependance of multistrange baryons. In peripheral collisions the yield of 
strange particles is small such that also here the canonical description should be applied. The canonical suppression of 
thermal particle phase-space increases with strangeness content of the particle. The exact conservation of strangeness 
requires that each particle carrying strangeness s has to appear e.g. with s other particles of strangeness one to satisfy 
strangeness neutrality condition. 



ho' 





FIG. 4. Particle multiplicities per participant normalized to its value in p+p system as a function of the number of 
participants A par t. calculated in the statistical model in the canonical ensemble. 

In Fig. 4 we calculate the multiplicity /participant of fi,S, and A relative to its value in a small system with 
only two participants |To[ . Thermal parameters were assumed here to be A part independent. Fig. 4 shows that the 
statistical model in (C) ensemble reproduces the basic features of WA97 data |l6|: the enhancement pattern and 
enhancement saturation for large A part indicating here that (GC) limit is reached. Fig. 4 also demonstrate different 
Apart dependence of strange and multistrange baryons. For small A par t this dependence is power like as describe by 
Eq.(22). The quantitative comparison of the model with the experimental data would require an additional assumption 
on the variation of [1b with centrality to account for larger value of B/B ratios in p+A than in Pb+Pb collisions 
Jl0|,[l6]]. The most recent results of NA57 showing an abrupt change of the enhancement for 3 are, however, very 
unlikely to be reproducible in terms of the canonical approach. 



V. SUMMARY AND CONCLUSIONS 



We have discussed the importance of the conservation laws in the application of the statistical model to the 
description of strangeness production in heavy ion collisions. We have presented the arguments that the more general 
treatment of strangeness conservation based on the canonical ensemble is required if one compares the model with 
experimental data for particle yields obtained in central A-A collisions at SIS energies or peripheral collisions at the 
SPS. In both situations the number of produced strange particles per event is still too small to use the asymptotic 
grand canonical ensemble. The time evolution of strangeness production and the approach to chemical equilibrium 
limit was discussed in the context of a kinetic approach. We have shown on few examples that the statistical model 
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predictions are consistent with the experimental data. A more complete presentation of the model versus data can 
be found in @^|||~t). 
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